Initial Conditions and Entanglement Sudden Death 
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We report results bearing on the behavior of non-local decoherence and its potential for being 
managed or even controlled. The decoherence process known as entanglement sudden death (ESD) 
can drive prepared entanglement to zero at the same time that local coherences and fidelity remain 
non-zero. For a generic ESD-susceptible Bell superposition state, we provide rules restricting the 
occurrence and timing of ESD, amounting to management tools over a continuous variation of initial 
conditions. These depend on only three parameters: initial purity, entanglement and excitation. 
Knowledge or control of initial phases is not needed. 



INTRODUCTION 



A continuously interacting background environment 
tends to destroy quantum state coherence. Quantum 
non-separability (entanglement), which plays a key role 
in various quantum computing and quantum communi- 
cation proposals [l|, is the best-known example of non- 
local coherence. Regarding separability, it is known that 
there is set of finite measure of separable states. If the 
steady state solution of a master equation for an initially 
entangled state lies on the interior of this set, then the 
transition from an entangled to a separable state has to 
occur after only a finite time of evolution. A version of 
these facts is discussed in specific terms in connection 
with an event which is called ESD (early-stage decoher- 
ence or more commonly entanglement sudden death) by 
Yu and Eberly Q and by Al Qasimi and James The 
inevitability of ESD for a class of initial states does not, 
however, make available the actual time when ESD will 
occur. Knowledge of the ESD time allows one, in at 
least some situations, to undertake preparations to avoid 
or delay its onset. In the following we address this ques- 
tion. 

Studies have been made to examine the ways in which 
ESD depends on dissipation mechanisms (e.g., ampli- 
tude damping [4], phase damping [5], and others l6[), 
and/or on interaction structures (e.g., interacting [7[ or 
non- interacting Q entangled partners, common (9[ or in- 
dependent [ToJ reservoirs), etc. Most of these factors are 
difficult to control in practice. In this Letter we extend 
the study of ESD control in a straightforward way by sys- 
tematically deriving connections between ESD and initial 
conditions, about which only isolated facts are noted, 
e.g. , in the study of ESD for Werner-like [ll| and X-type 
jl2| of initial states. We focus on the consequences of 
the initial preparation of entanglement, purity, and dou- 
ble excitation probability. Their initial values may be 
critical for an experiment's design and also for avoiding 
noise-induced ESD during the experimental time inter- 
val, but these initial values are not always easy to es- 
tablish precisely. However, for the generic case of ex- 
posure to amplitude noise, one can identify initial- value 
phase spaces which contain finite zones of safe prepa- 



ration. That is, within the well-defined boundaries of 
these zones the state evolution is guaranteed not to be 
ESD-susceptible. To our knowledge, this is the first indi- 
cation that such phase boundaries exist for any specific 
dissipation mechanism, and we provide formulas for the 
boundary curves. 

Additionally, analysis of the ESD-susceptible phases 
allows us to determine another feature of ESD that has 
remained generically unknown, namely the precise time 
when ESD will occur, if it is to occur at all. The ESD on- 
set time is found to be very sensitive to the initial purity 
in the vicinity of critical purity points. Especially for ini- 
tial states with very high or low entanglement, we show 
that a small amount of state impurity will accelerate the 
ESD process dramatically. Among potential advantages 
in practice, the conditions we have found allow protec- 
tion from ESD and may be useful in guiding experiments 
toward preparation of states that are immune to ESD or 
have relatively longer ESD times. 



II. PREPARED INITIAL STATE 

As an example we can consider the not quite per- 
fect preparation of Bell states |<£ + ) and |$ - ), which are 
known to be vulnerable to ESD in any combination. We 
suppose that a quantum computation task aims to use a 
generic combination such as 



\$ab) = cos(9|e)A|e)B + sin 0\g) A \g)B, 



(1) 



where \c)a and \q)a are orthogonal states of qubit A, etc., 
and may be interpreted as excited and ground states. In 
this state the degree of entanglement between the two 
qubits, as measured by concurrence [13j, is C = sin 26. 

In reality, in preparation, a state will deviate from the 
ideal target state. To sketch how this might occur we 
can begin with the state after its preparation. It is cou- 
pled to an environment that will cause post-preparation 
decoherence, and also still entangled with marginal en- 
tities denoted by M, i.e., forces, fields and objects that 
may participate in the preparation phase but then cease 
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interaction (see [3, EH)- Thus we write: 

|$(0)> = [cos0|e>A|e>B|mi) 

+ sin 6\g} A \g} B \m 2 } <8> \fa) a \fa)b- (2) 

Here # determines the degree of excitation of the two- 
party state via cos 2 — pn, sin 2 — ^44, and |0o)a an d 
|0o)b are the normalized initial states (usually ground 
states) of environmental reservoirs a and 6 respectively, 
and I mi) and 1 777,2) are normalized states of the marginal 
system M. Imperfect control of the preparation leads to 
a possibly mixed rather than pure initial state via the 
relation sin# cos (777,2 1 mi) = Pi4 < y/Pnpu- Here pij 
are the matrix elements of the initial two-qubit reduced 
density matrix 



/Wo) = 



/ pu pu \ 





\ p41 p 4 4 / 



(3) 




III. ENTANGLEMENT EVOLUTION 



FIG. 1: Time-dependent surfaces showing entanglement evo- 
lution parameterized in two ways: (left) C vs.r for a range 
of prepared double excitation values, and fixed purity; and 
(right) C vs.r for a range of purities for fixed double exci- 
tation. The symbols are defined in the text. In both plots 
the solid-color base-plane shows the region where C = 0, i.e., 
where ESD has already occurred. It is clear that ESD-free 
regions are accessible, but the extreme time-sensitivity of the 
ESD boundaries is not obvious here. 



the time dependent marginal states preserve the overlap 
relation (777,2 (t)\ m\ (t)) = (777,2 \mi). Then the two-qubit 
time-dependent reduced density matrix is obtained by 
tracing off the reservoir and marginal states, yielding: 



We next consider the dynamics of the two-qubit en- 
tanglement to see how it is affected by the initial state 
specification, and expose each qubit to the same ampli- 
tude damping process 



\e)A\(f>o)a 
\g)A\(/>o)a 



^Afo|e) A \fa) c 

|0>A|0O>a, 



/p^\9)a\(I>i)c 



(4) 



where \fa) a is another environmental reservoir state with 
a(fa\fa)a — 0, and q a is the probability that qubit 
A remains in its excited state. For the simplest illus- 
tration we assume exponential decay and take q a (t) = 
exp(— Ft) = exp(— r). Then the excitation-transfer prob- 
ability p a = 1 — q a will grow from to 1 irreversibly with 
increasing dimensionless time r. Qubit B and reservoir 
b are characterized similarly as (j4]) with decay probabil- 
ity q . Of course this is more specialized than needed 
- the rates of change Y need not be the same for both 
qubit s, and the amplitude channel can also be used to 
characterize non-dissipative evolutions such as produced 
by the XY spin interaction [16| or the Jaynes-Cummings 
interaction [17j. 

Amplitude damping takes the initial state (j2j) to the 
time-dependent state 



!*(*)> 



cos9^q a q h \ee)\^Q) a \<t>o)b\rni{t)) 
cos 0^/q a p b \eg)\^o) a \^i)b\mi(i)) 

COS0 v /Pa^|5 , e)|0i)a|0o)b|mi(t)) 

cos Oy/paPb I gg) | fa ) a \ fa ) b \ m x (t) ) 
sm6\gg)\fa) a \fa) h \m 2 {t)). (5) 



Here because of the fact that the marginal entity has 
ceased interaction with the qubit system after t = 0, 



Pab (t) = (6) 

/ PllQaQb Plly/qaQb 

pnq a Pb 

pllPaQb 

V P4i^/q a qb pllPaPb + p44 

The two-party entanglement (concurrence [l3|) obeys the 
X-state formula 1181: 



C(t) = 2^/q a (T)q b (T) 
x max{0, I P14 1 - 



Pn 



VPaCrRr)}, (7) 



which is graphed in two ways in Fig. [TJ 

Three parameters associated with the state mentioned 
above, that are presumably under good if not perfect 
preparation control [l9|, will now be selected to guide 
examination of evolution. These parameters are Co, pn 
and R, namely the initial entanglement (concurrence), 
initial excited state population, and a normalized value 
of purity. The true initial purity is given by 

P = l-2(pnp 4 4- |pi 4 | 2 ) 

= l-2 Pll (l-pii) + C 2 /2, (8) 

and the normalized version is defined as 

R = V2P- 1, (9) 

which takes values conveniently between and 1. 

IV. ESD SUSCEPTIBILITY PHASE DIAGRAMS 

The plots in Fig. [1] are only for a specific category 
of two-party states, but are systematically exhaustive in 
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FIG. 2: Entanglement phase diagrams of the initial two-qubit 
state in terms of the initial concurrence Co versus the double 
excitation number pn. The ESD phase (shaded area) and 
ESD-free phase (unshaded area) are separated by the dashed 
critical line p\\ — Co/2. 



FIG. 3: Entanglement phase diagram in terms of the ini- 
tial concurrence Co versus the initial normalized purity R. 
The ESD phase (shaded area) and ESD-free phase (un- 
shaded area) are separated by the dashed critical line R c = 
V2C 2 -2C + 1. 



the sense that any arbitrary combination of the Bell 
states is included, as well as any initial purity, excitation 
strength, and concurrence. It turns out that their ESD 
boundaries are directly analyzable, as follows. 

First we use (|6j) and ([7j) at r = 0, where q — 1 and 
p = 0, to characterize the entire physical domain in 
terms of pn and Co- Non- negativity of p requires the 
value of pn to be bounded between p™ ax and p™, i.e., 

1 ± yl — Cq /2, for any fixed Co. From the fact that 

Co < 1, one immediately has 



Pn S Co/ 2, < Pn 



(10) 



Within this domain, as r — >• oo, both p a and p^ and 
their product grow from to 1. Then it is easy to identify 
within the physical domain the boundaries between ESD 
and non-ESD phases. The condition C(r) =0 tells us 
that ESD occurs wherever 



y/PaPb > 



|P14| 
Pll 



Co 
2pn 



(11) 



From the fact that ^/p a Pb < 1, one immediately notes 
that ESD is restricted to the domain where the initial 
double excitation number is greater than half of the ini- 
tial concurrence, i.e., Co/2 < pn < p™ ax . The opposite 
region pff n < p n < C /2 is ESD-free. This leads to 
the phase diagram in Fig. [2j which locates the ESD and 
ESD-free phases in terms of initial concurrence Co and 
the double excitation probability p\\. The solid (in color, 
red and blue) boundary lines represent the values p™ ax 
and pif n as a function of Co- The dashed line pn = Co/ 2 
is the critical phase boundary where the state crosses 
from the ESD-free phase to the ESD-inevitable phase. 

An important implication of this analysis is that one 
needs to prepare the state with low double excitation, 
pn < Co/ 2 < 1/2, in order to avoid ESD. Therefore we 
will focus on this region of greatest interest and assume 
Pu < 1/2 in the following discussion. 



We now go a further step by employing the fact that 
the initial concurrence Co and the double excitation prob- 
ability pn are simply related to the initial state purity 
P as in (|8j). Then the ESD-onset condition ([TT]) can be 
rewritten in terms of the initial concurrence and normal- 
ized purity R as follows 



y/PaPb > 



Co 



VR 2 



r<2 1 



(12) 



where the non- negativity of the initial density matrix now 
requires that the normalized purity R obeys i? m in = Co 
and i^max = 1 for any fixed Co- Fig. [3] locates the ESD 
and ESD-free phases in Co-R space. The (in color, red 
and blue) solid lines represent the values i? max and R m in 
as functions of Cq. The dashed line is defined by 



Rc 



2Cl 



2C 



(13) 



which specifies the critical purity, as a function of Co, 
where the state crosses from the ESD-free phase to the 
ESD-inevitable phase. 

One sees that the states with very high or very low ini- 
tial concurrences are quite fragile and require very high 
purity in order to avoid ESD. States with initial concur- 
rence around 0.5 are the most robust, because they tol- 
erate the widest range 1] of variation of R or, in 
other words, have higher tolerance of preparation errors 
(permitting lower initial purity) before breaking down 
to be ESD-susceptible. This may have application in 
preparing ESD-free states. 



V. ESD ONSET TIME 

Although the phase boundaries are now located, there 
remains a question of near-boundary values of concur- 
rence. In the ESD-inevitable regions the onset of ESD 
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is not instantaneous so temporary preservation of entan- 
glement can be obtained if the parameter region close 
enough to the phase boundary for ESD onset can be well 
defined. 

Within an ESD-susceptible zone we denote by T the 
time when ESD occurs and entanglement permanently 
vanishes, i.e., the ESD onset time. That is, T is the 
value of r when relation ([T2]) is an equality. With p a = 



Pb = 1 



the onset time is then given by 



T = ln 



1 - Co - ^R 2 ~Cl 



(14) 



which is displayed in Fig. |4] as a function of both R 
and Co- One immediately sees from Fig. [4] that for any 
fixed initial concurrence Co G [0, 1], the ESD time is very 
sensitive to the initial purity in the vicinity of its critical 
values R c . Especially for the states having very high or 
very low initial concurrences, a small departure of the 
initial state from R = 1 will accelerate the ESD process 
dramatically. 




FIG. 5: ESD occurrence is identified here by time to onset of 
ESD. The graph can be read as a traffic signal, where green 
is a safe zone with no possibility of ESD, and red is the no-go 
zone where ESD is sure to occur quickly (within one life- 
time of the dissipative process or less), and the yellow zone 
is one where ESD is certain to occur but definitely deferred, 
i.e., ESD cannot occur before one lifetime, and so permits 
informed state preparation with the initial entanglement and 
purity values within its borders. 



VI. SUMMARY 
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FIG. 4: Entanglement sudden death time as a function of the 
initial entanglement Co and initial normalized purity R. The 
steepness of the boundary indicates existence of a narrow zone 
of "safe" state preparation, even though susceptible to ESD. 



An important practical point is that we can talk about 
three zones with clear meaning. One is the zone of infinite 
onset time, where ESD cannot occur, and one is the zone 
where ESD occurs inevitably and very rapidly. A more 
interesting zone is the one where the onset time for ESD is 
delayed in a definite way. A practical issue could be that 
we can tolerate one dissipative lifetime of decay, but only 
if we can be assured that ESD will not occur. This means 
that we need to know the tolerable equivalent parameter 
range for the initial preparation. To illustrate this case, 
we impose on Fig. [3] another boundary line, located at 
T = 1, and show the result in Fig. [5] The yellow zone 
shows how the parameter range available to the prepa- 
ration process can be expanded into the ESD-susceptible 
phase, tolerating some entanglement dissipation while re- 
maining certain to avoid ESD up to the predetermined 
time T = 1. 



In summary, we have reported new results bearing on 
the behavior of non-local decoherence and its potential 
for being managed or even controlled. As is well known 
and demonstrated [Sil ESD exists, i.e., decoherence 
processes can drive prepared entanglement to zero at the 
same time that easily monitored local coherences and fi- 
delity remain non-zero. Its inverse process, sudden birth 
of entanglement [21[ , is equally interesting and is the sub- 
ject of a later report. 

Until now there have been no rules of thumb or intu- 
itive guides giving reliable information about the likely 
occurrence or non-occurrence of ESD even in the simplest 
instances of entanglement and under the simplest of deco- 
herence mechanisms. We have considered a Bell superpo- 
sition state as practically the simplest possible entangle- 
ment scenario, and have subjected it to amplitude damp- 
ing, a well-understood decoherence mechanism. Similar 
steps have been taken previously but only for specific ini- 
tial density matrices generated ad hoc for demonstration 
purposes. 

In contrast, our findings amount to a first step toward 
defining a set of management tools that permit the range 
and extent of entanglement decoherence to be bounded 
and controlled. These tools, i.e., the values of initial 
purity, entanglement and excitation, are based only on 
knowledge of parameters that are presumably under good 
control during any specific entangled state preparation. 
More importantly, the zone-control nature of these man- 
agement tools, i.e., the fact that ESD susceptibility is de- 
termined over separate zones of a continuous variation of 
initial conditions, allows a well-defined amount of margin 
for errors in practical imperfect preparations. Notably, 
knowledge or control of initial phases is not needed. 
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